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Very energetic astrophysical events are required 
to accelerate cosmic rays to above 10 18 eV. 
Gamma-ray bursts (GRBs) have been proposed 
as possible candidate sources^HH. In the GRB 
fireball model, cosmic ray acceleration should be 
accompanied by neutrinos produced in the de- 
cay of charged pions created in interactions be- 
tween the high-energy cosmic-ray protons and 
gamma-rays 4 . Previous searches for such neutri- 
nos found none, but the constraints were weak as 
the sensitivity was at best approximately equal 
to the predicted flux^ED. Here we report an upper 



limit on the flux of energetic neutrinos associated 
with gamma-ray bursts that is at least a factor 
of 3.7 below the predictions^^ 121 . This implies 
that GRBs are not the only sources of cosmic rays 
with energies > 10 18 eV or that the efficiency of 
neutrino production is much lower than has been 
predicted. 



Neutrinos from GRBs are produced in the decay of 
charged pions produced in interactions between high- 
energy protons and the intense gamma-ray background 
within the GRB fireball, for example in the A-resonance 
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process p + 7 



n + 7r + . When these pions decay 



via 7r + — > and /i + — > e + j^ e i/ M , they produce a flux 

of high-energy muon and electron neutrinos, coincident 
with the gamma rays , and peaking at energies of sev- 
eral hundred Te\™JJ. Such a flux should be detectable 
using k m 3 -sca le instruments like the IceCube neutrino 
telescop e! 12 * 13 * (Suppl. Fig. [lj. Due to maximal mixing 
between muon and tau neutrinos, neutrinos from pion 
decay in and around GRBs will arrive at Earth in an 
equal mixture of flavors. We focus here only on muons 
produced in Vn charged-current interactions. As the 
downgoing cosmic ray muon background presents chal- 
lenges for the identification of neutrino-induced muons, 
we achieve our highest sensitivity for upgoing (northern 
hemisphere) neutrinos. However, the tight constraint of 
spatial and temporal coincidence with a gamma-ray burst 
allows some sensitivity even in the southern sky. One of 
the two analyses presented here therefore includes south- 
ern hemisphere gamma-ray bursts during the 59-string 
IceCube run. 

The results presented here were obtained while Ice- 
Cube was under construction using the 40- and 59-string 
configurations of the detector, which took data from 
April 2008 to May 2009 and from May 2009 until May 
2010, respectively. During the 59-string data taking pe- 
riod, 190 GRBs were observed and reported via the GRB 
Coordinates Network^, with 105 in the northern sky. 
Of those GRBs, 9 were not included in our catalog due 
to detector downtime associated with construction and 
calibration. Two additional GRBs were included from 
test runs before the start of the official 59-string run. 
117 northern-sky GRBs were included from the 40-string 
periocPto compute the final combined result. GRB po- 
sitions were taken from the satellite with the smallest re- 
ported error, which is typically smaller than the IceCube 
resolution. The GRB gamma-emission start (T s ^ ar ^) and 
stop (T s t p) times were taken by finding the earliest and 
latest time reported for gamma emission. 

As in our previous study^, we conducted two analyses 
of the IceCube data. In a model-dependent search, we 
examine data during the period of gamma emission re- 
ported by any satellite for neutrinos with the energy spec- 
trum predicted from the gamma-ray spectra of individual 
GRBs 6 3. The model-independent analysis searches more 
generically for neutrinos on wider time scales, up to the 
limit of sensitivity to small numbers of events at ± 1 day, 
or with different spectra. Both analyses follow the meth- 
ods used in our previous work 7 , with the exception of 
slightly changed event selection and the addition of the 
southern hemisphere to the model-independent search. 
Due to the large background of down-going muons from 
the southern sky, the southern hemisphere analysis is 
sensitive mainly to higher energy events (Suppl. Fig. 
[3]). Systematic uncertainties from detector effects have 
been included in the reported limits from both analyses 
and were estimated by varying the simulated detector 
response and recomputing the limit, with the dominant 
factor the efficiency of the detector's optical sensors. 
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FIG. 1. Comparison of results to predictions based on ob- 
served gamma-ray spectra. Th e sum med flux predictions 
normalized to gamma-ray spectra ^ 9 * 15 ! are shown in dashe d 
lines; the cosmic ray normalized Waxman-Bahcall AmP^^l [ s 
also shown for reference. The pre dicted neutrino flux, when 
normalized to the gamma rays^, is proportional to the ra- 
tio of energy in protons to that in electrons, which are pre- 
sumed responsible for the gamma-ray emission (ep/e-, here 
the standard 10). The flux shown is slightly modified^ from 
the original calculation 9 . (j> v is the average neutrino flux at 
Earth, obtained by scaling the summed predictions from the 
bursts in our sample (F v ) by the global GRB rate (here 667 
bursts/yeai^). The first break in the neutrino spectrum is 
related to the break in the photon spectrum measured by the 
satellites, and the threshold for photopion production, while 
the second break corresponds to the onset of synchrotron 
losses of muons and pions. Not all of the parameters used 
in the neutrino spectrum calculation are measurable from ev- 
ery burst. In such cases, benchmark valued were used for the 
unmeasured parameters. Data shown here were taken from 
the result of the model-dependent analysis. 



In the 59-string portion of the model-dependent anal- 
ysis, no events were found to be both on-source and 
on time (within 10° of a GRB and between T s t ar t and 
T s t p). From the individual burst spectra^ with the 
ratio of energy in protons vs. electrons e p /e e — 10 [Ref. 
E], 8.4 signal events were predicted from the combined 
2-year dataset and a final upper limit (90% confidence) 
of 0.27 times the predicted flux can be set (Fig. [l). This 
corresponds to a 90% upper limit on e p /e e of 2.7, with 
other parameters held fixed, and includes a 6% system- 
atic uncertainty from detector effects. 

In the model-independent analysis, two candidate 
events were observed at low significance, one 30 sec- 
onds after GRB 091026A (Event 1) and another 14 hours 
before GRB 091230A (most theories predict neutrinos 
within a few minutes of the burst). Subsequent examina- 
tion showed they had both triggered several tanks in the 
IceTop surface air shower array, and are thus very likely 
muons from cosmic ray air showers. In Fig. [2] are shown 
limits from this analysis on the normalization of E~ 2 
muon neutrino fluxes at Earth as a function of the size 
of the time window |At|, the difference between the neu- 
trino arrival time and the first reported satellite trigger 
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FIG. 2. Upper limits on E~ 2 power-law muon neutrino fluxes. 
Limits were calculated using the Feldman-Cousins method^ 
from the results of the model-independent analysis. The left 
y-axis shows the total number of expected f M events while 
the right-hand vertical axis (F„) is the same as in Fig. [I] A 
time window of At implies observed events arriving between 
t seconds before the burst and t afterward. The variation of 
the upper limit with At reflects statistical fluctuations in the 
observed background rate, as well as the presence of individual 
events of varying quality. The event at 30 seconds (Event 1) 
is consistent with background and believed to be a cosmic-ray 
air shower. 



time. As a cross-check on both results, the limit from this 
analysis on the average individual burst spectra^ during 
the time window corresponding to the median duration 
of the bursts in the sample (28 seconds) was 0.24 times 
the predicted flux, within 10% of the model-dependent 
analysis. 

Assuming that the GRBs in our catalog are a rep- 
resentative sample of a total of 667 per yeai^, we can 
scale the emission from our catalog to the emission of 
all GRBs. The resulting limits can then be compared 
to the expected neutrino rates from models that assume 
that GR Bs are the main sources of ultra high energy cos- 
mic ray d 4 ! 8 ! 1 ^ , with sampling biases of the same order 
as model uncertainties in the flux prediction s] 18 ^ 19 !. Lim- 
its from the model-independent analysis on fluxes of this 
type are shown in Fig. [3j 

These limits exclude all tested models^Hini with their 
standard parameters and uncertainties on those parame- 
ters (Figs. [l][3]). The models arc different formulations of 
the same fireball phenomenology, producing neutrinos at 
proton-photon (737) interactions in internal shocks. The 
remaining parameter spaces available to each therefore 
have similar characteristics: either a low density of high- 
energy protons, below that required to explain the cosmic 
rays, or a low efficiency of neutrino production. 

In the GRB fireball, protons are believed to be acceler- 
ated stochastically in collisions of internal shocks in the 
expanding GRB. The neutrino flux is proportional to the 
rate of 757 interactions, and so to the proton content of the 



FIG. 3. Compatibility of cosmic-ray flux based models with 
observations. Shown are the allowed values of neutrino flux 
vs. the neutrino break energy (eb) in comparison to model pre- 
dictions with estimated uncertainties. Data were taken from 
the model-independent analysis from the time window corre- 
sponding to the median duration of the GRBs in our catalog 
(|At| = 28 seconds). Spectra are represented here as broken 
power laws (cf> v ■ {E^ 1 /eb, E < eb; E~ 2 , E > eb}) with a break 
energy eb corresponding to the A resonance for P7 interac- 
tions in the frame of the shock. The muon flux in IceCube is 
dominated by neutrinos with energies around the first break 
(eb). As such, the upper break, due to synchrotron losses of 
7r , has been neglected here, as its presence or absence does 
not contribute significantly to the muon flux and thus does 
not have a significant effect on the presented limits, eb is re- 
lated to the bulk Lorentz factor F (eb oc T 2 ); all of the models 
shown assume T ~ 300. The value of T corresponding to 10 7 
GeV is > 1000 for all models. Vertical axes are related to 
the accelerated proton flux by the model-dependent constant 
of proportionality / w . For models assuming a neutron-decay 
origin of cosmic rays (RacherP and Ahlers^) is indepen- 
dent of F; for others (Waxman-BahcalP) f n oc T~ 4 . Error 
bars on model predictions are approximate and were taken 
either from the original papers, where included^, or from the 
best-available source in the literature^ otherwise. The errors 
are due to uncertainties in f n and in fits to the cosmic-ray 
spectrum. Waxman-BahcalP and RacherP fluxes were calcu- 
lated using a cosmic ray density of 1.5 — 3 x 10 44 erg Mpc -3 
yr" 1 , with 3 x 10 44 the central valucP. 

burst by a model-dependent factor. Assuming a model- 
dependent proton ejection efficiency, the proton content 
can in turn be related to the measured flux of high-energy 
cosmic rays if GRBs are the cosmic ray sources. Limits 
on the neutrino flux for cosmic ray normalized models 
are shown in Fig. [3j each model prediction has been 
normalized to a value consistent with the observed ultra 
high-energy cosmic ray flux. The proton density can also 
be expressed as a fraction of the observed burst energy, 
directly limiting the average proton content of the bursts 
in our catalog (Fig. E|. 

An alternative is to reduce the neutrino production ef- 
ficiency, for exa mple b y modifying the physics included 
in the prediction s] 19 ! 20 ! or by increasing the bulk Lorentz 
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Bulk Lorentz Factor (D 



FIG. 4. Constraints on fireball parameters. The shaded re- 
gion, based on the result of the model-dependent analysis, 
shows the values of GRB energy in protons a nd t he average 
fireball bulk Lorentz factor for modeled fireballs^' allowed by 
this result at the 90% confidence level. The dotted line in- 
dicates the values of the parameters to which the completed 
IceCube detector is expected to be sensitive after 3 years of 
data. The standard values considered^ are shown as dashed- 
dotted lines and are excluded by this analysis. Note that the 
quantities shown here are model-dependent. 



boost factor T. Increasing T increases the proton en- 
ergy threshold for pion production in the observer frame, 
thereby reducing the neutrino flux due to the lower pro- 
ton density at higher energies. Astrophysical lower limits 
on r are established by pair production arguments^, but 
the upper limit is less clear. Although it is possible that 
r may take values of up to 1000 in some unusual bursts, 
the average value is likely lower (usually assumed to be 
around 300^1) and the non-thermal gamma-ray spectra 
from the bursts set a weak constraint that T < 2000^2. 
For all considered models, with uniform fixed proton con- 
tent, very high average values of T are required to be 
compatible with our limits (Figs. [3j [4j) . 

In the case of models wh ere c osmic rays escape from 
the GRB fireball as neutroniPE^, the neutrons and neu- 
trinos are created in the same pj interactions, directly 
relating the cosmic ray and neutrino fluxes and remov- 
ing many uncertainties in the flux calculation. In these 
models, T also sets the threshold energy for production 
of cosmic rays. The requirement that the extragalactic 
cosmic rays be produced in GRBs therefore does set a 
strong upper limit on T: increasing it beyond ~ 3000 
causes the proton flux from GRBs to disagree with the 
measured cosmic ray flux above 4 x 10 18 eV, where extra- 
galactic cosmic rays are believed to be dominant. Limits 
on r in neutron-origin models from this analysis (> 2000, 
Fig. |3| are very close to this point, and as a result all 
such models in which GRBs are responsible for the entire 
extragalactic cosmic-ray flux are now largely ruled out. 

Although the precise constraints are model dependent, 



the general conclusion is the same for all the versions of 
fireball phenomenology we have considered here: either 
the proton density in gamma ray burst fireballs is sub- 
stantially below the level required to explain the highest 
energy cosmic rays or the physics in gamma ray burst 
shocks is significantly different from that included in cur- 
rent models. In either case, our current theories of cos- 
mic ray and neutrino production in gamma ray bursts 
will have to be revisited. 
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IceCube Lab 




SUPPL. FIG. 1. The IceCube neutrino observatory. The Ice- 
Cube detector instruments a volume of 1 km 3 of glacial ice 
at the South Pole, sensitive to neutrinos of TeV and higher 
energjD^l (Suppl. Fig. [2]). Neutrinos are detected by observ- 
ing Cherenkov light emitted by secondary charged particles 
produced in neutrino-nucleon interactions- 1311 , and their arrival 
direction is obtained from the timing pattern of the detected 
light. The finished detector is composed of 5160 digital optical 
modules (DOMs), each containing a 10-inch photomultiplier, 
with 60 placed at depths between 1450 and 2450 m on each 
of 86 vertical strings. IceCube is complemented by a sur- 
face air shower array called IceTorP 1 , with two tanks located 
above each of the IceCube strings. The colors at the top in- 
dicate the detector at various stages of deployment. IceCube 
achieves its best angular resolution for muons produced in 
charged-current interactions (0.6° for E v > 100 TeV). Com- 
bined with the increased detector effective volume afforded 
by the long distances traveled by the secondary muons, such 
events usually provide the highest sensitivity for searches for 
neutrino point sources. 
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SUPPL. FIG. 2. Effective area of the IceCube neutrino tele- 
scope using the event selections of the model-dependent and 
model-independent analyses, averaging over the 40- and 59- 
string detector configurations and zenith angles according to 
the distribution of bursts in the catalog. The effective area 
of the model-independent event selection is in general some- 
what larger, due to using a weight scheme instead of hard 
cuts - however, the extra events so included are typically low 
quality and so have low weights when computing final re- 
sults. The model-independent average effective area includes 
the southern hemisphere for the 59-string portion of the anal- 
ysis (Suppl. Fig. [3]|. 

EFFECTIVE AREAS 

The detector effective areas (Suppl. Figs. can 
be used to estimate the detector response for an arbi- 
trary neutrino flux. Convolution of a flux with the ef- 
fective area will give the expected event rate in IceCube. 
Presented effective areas are the average of the effective 
areas for muon neutrinos and muon antineutrinos and 
correspond to the expectation value of the detector ef- 
fective area under variations to account for systematic 
uncertainties in the detector simulation. The increase 
in effective area between the 40- and 59-string detector 
configurations is due to the 50% increase in geometrical 
area of the detector, a more favorable detector geometry, 
and improvements in the event selection and reconstruc- 
tion techniques (Suppl. Fig. Eh. Data files containing 
all the effective areas plotted here are included in the 
supplementary information (Suppl. Tables 1-6). 



The GRB catalog used in this analysis was syn- 
thesized from GCN notices and can be obtained us- 
ing the GRBweb database available at |http : / /grbweb . 
icecube.wisc.edu/. IceCube-40 operated from April 
5, 2008 until May 20, 2009 and IceCube-59 operated 
from May 21, 2009 until May 31, 2010. GRB090422 
and GRB090423, though before the official 59-string start 
date, occurred during test runs of the 59-string detector 
and so are included in the 59-string catalog. 



COMBINATION OF DATASETS 

The results presented use a combination of the IceCube 
40- and 59-string datasets. In both analyses, all GRBs 
were individually simulated and this simulation was ap- 
plied to the detector running at the time of the GRB. The 
simulated events from the full GRB catalog were treated 
as a combined dataset, which was then compared to the 
combined result from both detector configurations. Sys- 
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SUPPL. FIG. 3. Average effective area of the IceCube neu- 
trino telescope in its 59-string configuration to fluxes of muon 
neutrinos from the northern and southern hemispheres, us- 
ing the event selection from the model-independent analysis. 
At low energies, the ability to use the Earth to filter out 
neutrinos from cosmic ray air showers reduces backgrounds 
and improves the sensitivity of the detector to neutrinos from 
the northern sky. As the neutrino energy increases, so does 
the neutrino-nucleon cross-section, increasing the neutrino in- 
teraction probability and the chances of a detection in Ice- 
Cube. At very high energies, backgrounds from cosmic-ray 
muons are substantially reduced and the neutrino-nucleon 
cross-section becomes large enough that neutrinos from the 
northern hemisphere will be absorbed in the Earth before 
reaching the detector, making the southern hemisphere the 
region of highest sensitivity for E v > 1 PeV. Note that the 
effective area is not a direct estimate of the sensitivity, due 
to variable backgrounds as a function of zenith angle and en- 
ergy. In general, backgrounds are highest in the southern 
hemisphere and at low energies. 



SUPPL. FIG. 4. Comparison of IceCube effective area be- 
tween the 40- and 59-string detectors in the northern hemi- 
sphere using the model-independent event selection. The in- 
crease in effective area is related to the increased geometric 
area of the detector, as well as to improved event selection 
and reconstruction techniques. Low-energy sensitivity was 
especially enhanced as a result of a more favorable detector 
geometry and the deployment of the first of the strings of the 
Deep Core subdetector. As a result, the integrated effective 
area of the 59-string detector for this analysis is more than 
1.5 times that of the 40-string detector. 



tcmatic uncertainty estimates computed by variation of 
detector parameters were likewise applied separately to 
each detector configuration and GRB sub-catalog, and 
then used in a combined fit for the determination of the 
final limits. 



